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Abstract. We investigate the cold nuclear matter effects on J/tp production, whose 
understanding is fundamental to study the quark-gluon plasma. Two of these effects 
are of particular relevance: the shadowing of the parton distributions and the nuclear 
H absorption of the cc pair. If J/-0's are not produced via a 2 — > 1 process as suggested 

by recent theoretical works, one has to modify accordingly the way to compute the 
nuclear shadowing. This naturally induces differences in the absorption cross-section 
fit to the data. A careful analysis of these differences however requires taking into 
account the experimental uncertainties and their correlations, as done in this work for 
dAu collisions at ^/satjv = 200 GeV, using several shadowing parametrisations. 
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1. Introduction 

Relativistic nucleus-nucleus {AB) collisions are expected to produce a deconfined state 
of QCD matter - the Quark Gluon Plasma (QGP) - at high enough densities or 
temperatures. It has long been suggested [1] that the J/ip meson would be sensitive 
to Hot and Dense Matter (HDM) effects, through mechanisms like the dissociation of 
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the cc pair due to the colour Debye screening. A significant suppression of the J /if) 
yield was observed by the PHENIX experiment in CuCu [2] and AuAu [3] collisions 
at a/stvjv = 200 GeV. However, before giving any interpretation, Cold Nuclear 
Matter (CNM) effects have to be properly disentangled and subtracted. They are 
known to impact the J ft/) production in proton-nucleus (pA) or deuteron-nucleus (dA) 
collisions, where the deconfinement conditions can not be reached. Such non-trivial 
effects are demonstrated by the PHENIX dAu data [4] obtained at the same energy. 
Two CNM effects are of particular importance [5]: (i) the shadowing of the initial 
parton distributions (PDFs) due to the nuclear environment, and (ii) the breakup of 
cc pairs consecutive to multiple scatterings with the remnants of the projectile and target 
nuclei, referred to as the nuclear absorption. Recent theoretical works incorporating 
QCD corrections or s-channel cut contributions have emphasized [6, 7] that the Colour- 
Singlet (CS) mediated contributions are sufficient to describe the experimental data for 
hadroproduction of both charmonium and bottomonium systems without the need of 
Colour-Octet (CO) contributions. Furthermore, recent works [8] focusing on production 
at e + e~ colliders have posed stringent constraints on the size of CO contributions, 
which are the precise ones supporting a 2 — > 1 hadroproduction mechanism [9]. As a 
consequence, J /if) production at low and mid Pt likely proceeds via a 2 — > 2 process, 
such as g + g — > J /if) + g, instead of a 2 — > 1 process. As we have shown in previous 
studies [10, 11, 12], this modifies both the way to compute the nuclear shadowing and 
its expected impact on the J / if) production. In this work, we shall focus on the changes 
induced on the rapidity dependence of the J/ if) nuclear modification factor in <iAu, while 
using several parametrisations of the nuclear PDF. As in [4] - where CNM effects were 
computed based on a 2 — > 1 kinematics - we shall use the same <iAu data to derive 
the absorption cross-section <7 a bs required on top of the shadowing, but assuming here a 
2—7-2 underlying partonic process. We shall compare the results found in both schemes, 
with a special emphasis on the limitations from the experimental uncertainties. 

The article is organized as follows. In section 2, we will describe our model and the 
method chosen to carry a data-driven evaluation of the nuclear absorption cross-section. 
And in section 3, we will present and discuss our results before concluding. 

2. Our approach 

To describe the J /if) production in nucleus collisions, our Monte Carlo framework [10, 13] 
is based on the probabilistic Glauber model, the nuclear density profiles being defined 
with the Woods-Saxon parameterisation for any nucleus A > 2 and the Hulthen 
wavefunction for the deuteron [14]. The nucleon-nucleon inelastic cross section at 
t/snn = 200 GeV is taken to ctnn = 42 mb and the maximum nucleon density to 
p = 0.17 nucleons/fm 3 . 

In order to study the J / if) production, we need to implement in our Monte Carlo the 
following ingredients: the partonic process for the cc production and the CNM effects. 
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2.1. Partonic process for the cc production 

Most studies on the J/ip production in hadronic collisions rely on the assumption that 
the cc pair is produced by the fusion of two gluons carrying some intrinsic transverse 
momentum kx- The partonic process being a 2 — > 1 scattering, the sum of the gluon 
intrinsic transverse momentum is transferred to the cc pair, thus to the J/if> since 
the soft hadronisation process does not alter significantly the kinematics. This is 
supported by the picture of the Colour Evaporation Model (CEM) at LO (see [9] and 
references therein) or of the CO mechanism at a 2 [15]. In such approaches, the transverse 
momentum Pt of the J /if} entirely comes from the intrinsic transverse momentum of 
the initial gluons. 

However, the average value of k T is not expected to go much beyond ~ 1 GeV. 
So this process is not sufficient to describe the Pt spectrum of quarkonia produced in 
hadron collisions [9] . For Pt > 2 — 3 GeV, most of the transverse momentum should 
have an extrinsic origin, i.e. the J/^'s Pt would be balanced by the emission of a 
recoiling particle in the final state. The J/if) would then be produced by gluon fusion 
in a 2 2 process with emission of a hard final-state gluon. This emission, which is 
anyhow mandatory to conserve C-parity, has a definite influence on the kinematics of 
the J /if) production. Indeed, for a given J/if> momentum (thus for fixed rapidity y and 
Pt), the processes discussed above, i.e. g + g — » cc — >■ J /if) (+X) and g + g — >■ J /if} + g, 
will proceed on the average from initial gluons with different Bjorken-x. Therefore, they 
will be affected by different shadowing corrections. From now on, we will refer to the 
former scenario as the intrinsic scheme, and to the latter as the extrinsic scheme. 

In the intrinsic scheme, we use the fits to the y and Pt spectra measured by 
PHENIX [17] in pp collisions at ^snn = 200 GeV as inputs of the Monte-Carlo. Indeed, 
the measurement of the J/if> momentum completely fixes the longitudinal momentum 
fraction carried by the initial partons: 

xi, 2 = -^=exp (±y) = x° 12 (y,P T ), (1) 

y/SNN 

with the transverse mass tut = ^JM 2 + P 2 , M being the J/ip mass. 

On the other hand, in the extrinsic scheme, information from the data alone - the 
y and Pt spectra - is not sufficient to determine x± and x 2 . Actually, the presence 
of a final-state gluon introduces further degrees of freedom in the kinematics, allowing 
several (£1,2:2) for a given set (y,Pr). The four- momentum conservation explicitely 
results in a more complex expression of x 2 as a function of (xi, y, Pt)'- 

Xim T ^s NN e- y - M 2 

2 y/S N N(y/S N NXl ~ m T e V ) 

Equivalently, a similar expression can be written for function of (x 2 , y, Pt)- Even 

if the kinematics determines the physical phase space, models are anyhow mandatory 
to compute the proper weighting of each kinematically allowed (x±,X2)- This weight 
is simply the differential cross section at the partonic level times the gluon PDFs, 
i.e. g(xi, HF)g(x2, Hf) do-gg->j/ip+g/dy dPx dxidx2- In the present implementation of our 
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code, we are able to use the partonic differential cross section computed from any 
theoretical approach. For now, we use the one from [6] which takes into account the 
s-channel cut contributions [16] to the basic CS model and satisfactorily describes the 
PHENIX pp data [17] down to very low P T [10]. 

2.2. Shadowing and nuclear absorption 

To obtain the J / ip yield in pA and AA collisions, a shadowing-correction factor has to 
be applied to the J/ip yield obtained from the simple superposition of the equivalent 
number of pp collisions. This shadowing factor can be expressed in terms of the ratios 
Rf of the nuclear Parton Distribution Functions (nPDF) in a nucleon belonging to a 
nucleus A to the PDF in the free nucleon: 



The numerical parameterisation of Rf(x,Q 2 ) is given for all parton flavours. Here, 
we restrict our study to gluons since, at high energy, the J/ip is essentially produced 
through gluon fusion [9] . In order to see how the CNM effects can vary depending on the 
different shadowing parametrisations used as an input, we will consider three of them: 
nDSg [18] at LO, EKS98 [19] and EPS08 [20]. They span the current evaluation of the 
uncertainty [21] on the gluon nPDF, from a small to a very large antishadowing. 

The second CNM effect that we are going to take into account concerns the nuclear 
absorption. In the framework of the probabilistic Glauber model, this effect is usually 
parametrised by introducing an effective absorption cross section cr a b s of the pre-resonant 
cc pair when propagating in the nuclear medium. In the following, we shall compare the 
data-constrained values of cr a b s given two different partonic cc production mechanisms 
- intrinsic and extrinsic - and three shadowing parametrisations as cited above. 

2.3. A data-driven evaluation of the nuclear absorption cross-section 

We present here the derivation of the <T a b s values consistent with the PHENIX dAu data, 
taking into account the various experimental uncertainties and their correlations. To 
do so, we have chosen to follow the procedure from [4, 22]. For a given choice of the cc 
production mechanism and nPDF parametrisation, the best possible agreement of the 
theory to the data is obtained for the value of <r a b s that minimizes the quantity: 



with (Tj = (jj [di + CfoO"^ + e c dicr c ), di being the set of experimental values from PHENIX, 
/i, the respective values predicted by the theory for a given set of parameters p (i.e. 
the nPDF parametrisation and cr a b s ), cr« are the point-to-point uncorrelated errors 
(statistical and systematic), the point-to-point correlated systematic errors, a c the 
global systematic error on the normalisation of the data and e^ ;C ) the fractions of the 
systematic uncertainties cr( 6 . c ) used to shift the data points d^. 




(3) 
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3. Results and discussion 

In the following, we present our results for the J/ip nuclear modification factor in dAu 
collisions: -R^Au = dN^j^/ (N co \\)dNpl^ , where dN^^(dNpl^) is the observed J/ip yield 
in dAu (pp) collisions and (A co n) is the average number of nucleon-nucleon collisions 
occurring in one dAu collision. Without nuclear effects, -R^Au should equal unity. 
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Table 1: cr a b s extracted from fits of -RdAu vs y for the intrinsic (left) and extrinsic (right) 
schemes, when considering all the different types of errors on the data, together with 
the corresponding x 2 obtained for the best fit. 

In Table. 1, we recall the results from [12], with the value of cr a b s corresponding to 
the best fit to PHENIX data RdAu vs y, and the one standard deviation uncertainties. 
As in [4], this extraction relies on the assumption that <r a b s is independent of y. The best 
agreement to the data is obtained in the extrinsic scheme with EPS08. The larger is 
the antishadowing, the larger are the differences between both schemes. This is visible 
in the \ 2 f° r the best fit or on the different shapes of the curves in Fig. 1 which shows 
the obtained CNM effects together with PHENIX data for -R^Au vs y. In this plot, the 
anti-shadowing peak in -RdAu is systematically shifted towards larger y for the extrinsic 
scheme with respect to the one in the intrinsic case. This reflects the larger value of 
the gluon momentum fraction %2 in the Au nucleus needed to produce a J/ip when the 
momentum of the final state gluon is indeed accounted for. We shall now focus on the 
bands picturing the CNM effects in both schemes which account for the experimental 
uncertainties. For the time being, the current size of the errors unfortunately does not 
allow to distinguish between the two approaches if assuming a constant cr a b s with y. 

Using the same procedure, we have qualitatively studied three scenarios with thirty 
times more statistics and a) no improvement, b) a reduction of 35% and c) a reduction of 
50% of the systematics. We have found that, with the same assumption on cr abs (namely 
constant vs y), the latest dAu data (2008) will not be sufficient to distinguish between 
a 2 1 and a 2 — > 2 production mechanism by only looking at the ^-dependence of 
RdAu, unless the (anti) shadowing is as strong as encoded in the EPS08 parametrisation. 
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(a) nDSg at LO. (b) EKS98. (c) EPS08. 



Figure 1: RdAu vs y for the intrinsic and extrinsic schemes and for the three nPDF 
parametrisations, compared to PHENIX data [4]. For each scheme, the central band 
represents the range in cr a b s consistent with the data within one standard deviation, 
when taking into account the point-to-point uncorrelated (bar) and correlated (box) 
errors. The corresponding <j a b s values are reported in the legend. The outer band is the 
obtained range in <7 a b s when the global error on the data normalisation is also considered. 
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